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Abstract. We obtain exact results for the effective diffusion constant of a two 
dimensional Langevin tracer particle in the force field generated by charged point 
scatterers with quenched positions. We show that if the point scatterers have a screened 
Coulomb (Yukawa) potential and are uniformly and independently distributed then the 
effective diffusion constant obeys the Volgel-Fulcher-Tammann law where it vanishes. 
Exact results are also obtained for pure Coulomb scatterers frozen in an equilibrium 
configuration of the same temperature as that of the tracer. 
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1. Introduction 

The bulk transport properties of random media are of great importance in physics and 
engineering. The computation of the bulk diffusivity, conductivity, permeability and 
dielectric constant from the statistical properties of their local fluctuations is referred to 
as homogenisation. These bulk quantities describe the transport properties over large 
length and time scales and will normally have well defined values when the statistics of 
the local fluctuations are homogeneous under translation. We will consider situations 
that are also statistically isotropic since a breakdown of isotropy involves further levels 
of complexity, important and interesting though they are. 

In this paper we compute the effective diffusion constant for tracer particles 
travelling through a medium of fixed but randomly distributed centres of force. We 
we assume for simplicity that the centres of force are all identical but may have either 
positive or negative charge with equal probability. The tracer therefore experiences a 
randomly fluctuating potential with zero mean. This model of a disordered medium is 
a reasonable, if simplified, description of many real physical systems where the random 
potential is induced by impurities inserted in a regular background. 

Many studies of diffusion in a random potential have been based on the assumption 
that it can be represented by a Gaussian field [TJ El El HI E] Diffusion in a non- 
Gaussian field has been studied using perturbation theory 0. However the calculation 
is significantly more complicated than that for the Gaussian case. The significance of 



Effective diffusion constant in a two dimensional medium of charged point scatterers 2 



our investigation is that we obtain exact results for a class of problems, diffusion in 
random potentials of point scatterers, that are non-Gaussian and for which few results 
are known. The results are exact in one and two-dimensions but experience suggests 
that they may be indicative for higher dimensional cases. 



2. The Model 

The position X i: of a particle, of negligible inertia, subject to a white noise, Wi(t), and 
a potential 4>(x), satisfies the equation 

X i = w i {t) + X Vi<f>{X) , (1) 

where 

(w i (t)w j (lf)) = 2K 6 ij 6(t-lf) , (2) 

and (• • •) denotes an average over the white noise. The Einstein relation implies that 
the local, or bare, diffusivity k and the coupling to the potential gradient A are related 
by the equation 

where T is the absolute temperature in appropriate units. The probability density, 
p(x,t), for the position of the particle obeys 

= Ko V 2 p + A V • (pVtf) . (4) 

at 

The effective diffusivity, appropriate to the dispersion of the particle at large times 
and distances is 

K (9) = J_ lim (231 (5) 

2£>*-oo t K J 

where D is the dimension of space. The mean squared displacement of the particle is 
given by 

= jd D xx 2 P { X) t) . (6) 

The superscript g is used to denote the fact that we are considering the effective diffusion 
constant for a particle diffusing in a gradient field. 

The description of the model is completed by specifying the structure of the 
potential (j>(x) in terms of the scattering centres and their charges. In this paper we will 
consider potentials of the type 

N 

^( X ) = Yl <lnV{x - X n ) , (7) 

71=1 

generated by iV scatterers frozen in the volume V. Here q n — ±1 is the charge on the 
scatterer n, x n is its position and V(x) the potential at the point x due to a positively 
charged background particle at the origin. In this paper we shall consider distributions of 
the x n and q n such that the field is statistically invariant under global transformation 
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— > —cf). This is clearly the case for distributions of the x n which are homogeneous 
and isotropic along with the condition that the distribution of the q n is invariant under 
the global charge transformation q n — > —q n . An example of this type of distribution is 
one where the positions x n are distributed uniformly and independently in the volume 
V and the charges are taken to be independent and ±1 with probability \. In this case 
the disorder average is given by 

(0) d = U E Wj^n O . (8) 

One could also consider an strictly electroneutral ensemble where N/2 of the scatterers 
have the charge ±1, here the disorder average is thus 

(°)<i = II \ I d ° X n O , (9) 



n 



we shall see however in the thermodynamic limit these two problems have the same 
diffusion constant. Another realisation of an ensemble of the type mentioned above is 
one where a system of N/2 positive and negative charges, interacting via the Hamiltonian 

N 

H = E QiQjV{xi - xj) , (10) 

i<j 

are allowed to interact and equilibriate at some inverse temperature /3', they are then 
frozen in this equilibrium configuration giving a distribution of the x n of 

N 



P({x, l }) = ^ r) ex V ^-P'Y, m V(x t -x J )) , (11) 



where Z(fi') is the partition function for the system. Clearly the uniform distribution 
of equation © is recovered in the limit /?' — > 0. 



3. Associated permeability model 

In a previous paper we showed that there is a strong connection between our gradient 
flow model and the problem of computing the effective long range diffusivity, of 
particle with a locally random diffusivity field n(x) . A Langevin particle diffusing in a 
random medium with local diffusivity n(x) has a probability density which satisfies 

= V ■ k(x)V P . (12) 

In this case also the long time transport can also be partially described by the 
asymptotic behaviour of the mean squared displacement. 

(X 2 t ) ~ 2D^t , (13) 

where reW is the bulk or effective diffusivity of the medium. 

The problem of calculating /4 p) is also equivalent calculating the bulk permeability 
of a porous medium, where the flow is described by Darcy's law [7j , given that the local 
permeability is k(x). It is also equivalent to calculating the effective conductivity or 
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effective dielectric constant in media where k(x) represents the local conductivity or 
dielectric constant. The superscript p here is thus used to denote that the effective 
diffusion constant corresponds to an effective permeability problem. The problem of 
calculating Kg' has been considered by numerous authors via exact relations in one 
and two dimensions [3 El El EH] and via perturbation techniques in three dimensions 

[m na nu 

The relevant connection between the permeability problem and the gradient flow 
problem emerges when the local diffusivity or permeability is chosen to be 

k(x) = k exp ((3<f>(x)) . (14) 



In [l2j it was shown that the bulk permeability or effective diffusion constant for equation 
(fT2J) with a local diffusivity given by equation ()14j) and that for the gradient flow problem 
equation (j3J with the same field cf> are related by 

4 9) «o (15) 



(p) 7g 



where 



k = — J d D x k(x) , (16) 

The relation equation (|15[) holds in all dimensions where the corresponding n e exist. In 
a statistically homogeneous system with short range correlations in the field 0, in the 
limit of large V we expect k to be self averaging, or realisation independent, i.e. 

1 r*. 



K = («) d = (— J d x K(x)) d , (17) 
4. Duality relation and the two dimensional result 

For two-dimensional systems we can obtain exact results by combining equation (|15|) 
with a standard duality theorem [HJ El QUI- For completeness we derive the relevant 
result in a form appropriate to the problem under consideration. 

As a preliminary we note that when the fields <j>{x) and —(f)(x) are statistically 
equivalent it follows that local diffusivities k(x) and k'(x) = Kq/k(x) are also 
statistically equivalent. We view equation (fT2|) as describing the flow of the density 
p(x) through the random medium. The associated current is j(x) = n{x)'S7p{x) . In the 
steady state we have 

V • (k(z)Vp(z)) = , (18) 

and on averaging over the ensemble of samples we have, assuming a constant mean 
current and density gradient, 

(j( x )) = (k(x)Vp(x)} = K^iVpix)} . (19) 

In two-dimensions equation (|18J) implies that there exists a dual field x( x ) such that 

j(x) = n(x)Vp(x) = K n x Vx(^) , (20) 
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where n is a unit vector orthogonal to the two-dimensional plane of the problem. We 
also have the dual equation 

k'(x)Vx(x) = —K Q n x Vp(x) (21) 

which implies 

V-(k'(x)V X (x)) = . (22) 

Since k'{x) is statistically equivalent to k{x) it follows that 

(k\x)Vx(x)) = ^(Vx(x)) , (23) 

with the same effective permeability , as in the original problem. We have then the 
two dual results 

4 p) (Vp(x)) = n Q n x (V*(x)) , (24) 

and 

k { p\Vx{x)) = -KqU x (Vp(x)) . (25) 
The consistency of these equations implies that 

4 p) = «o • (26) 
Combining this result with equation (fTB^) we obtain the general result valid in two 
dimensions 

«?> = i , (27) 

K 

This equation lj27jl was previously used by the authors to solve the problem of calculating 
in a homogeneous Gaussian field in two dimensions to give 

«W = Koexp(-^A(0)) , (28) 

where A is the connected two-point correlation function defined by 

A(x - x') = (<P{x)<P{x')) d - (<P{x)) d (<P{x')) d , (29) 

the average again being over the sample disorder. This result is particularly interesting 
as it showed that the renormalisation group result for PQ|2] in dimension D, given 
by 

#(i2G?) = « expf-^A(0)), (30) 



D 

was exact in two dimensions. The equation (|3U|) was already known to be exact in one 
dimension [2]. Indeed in one dimension the effective permeability is given for a general 
homogenous k(x) by 

# =(!)"' (3!) 

i.e. the harmonic mean of the local permeability. Using the result equation (|15|) we 
thus obtain the solution of the corresponding gradient flow problem in one dimension 
to be 

Ki 9) . , / T 



e 



(32) 
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5. Some Examples 

For the problem of tracer particles moving through scatterers uniformly and 
independently distributed in two dimensions, we can combine equation (jHJ) with equation 
(fTTjl to obtain 



n 



jd^n^j2exp((3 qi V(-x n )) = \l + ±Jd 2 x (cosh(/3V(a;))-l 

<?n J 



(33) 



We now take the thermodynamic limit, that is N — > oo at fixed p = N/V and obtain 
the exact result 



K exp (-p J d 2 x (cosh(/?V(x)) - 1)) . (34) 

The distribution given by equation (JHJ) for an exactly neutral system of scatterers of 
density p also leads to the same result equation (J34j) in the thermodynamic limit. 

It is interesting to review the relationship of this result to equation (J2E1) for the 
model in which <p(x) is a Gaussian field PJ Ej with zero mean. In the case of uniformly 
distributed randomly charged scatterers the two-point correlation function of the field 
is given by 

A(x -x')=p [ d 2 x" V(x - x")V(x' - x") . (35) 



Applying the Gaussian "prediction" equation (j28j) for the diffusivity therefore gives 

= K exp j-^p J d 2 x V 2 (x)| . (36) 

Comparing equations (|3l)j) and (JHU) we see that the result for the Gaussian field is 
recovered from that of the point scatterers in the limit /3 — ^ while pj3 2 = c with c 
a constant. This limit where the scatterer disorder acquires a Gaussian character can 
also be found by examining the perturbation theory for this problem and holds in all 
dimensions. It is interesting to note that the Gaussian limit does not necessarily hold, 
as one may have naively expected from the central limit theorem, simply in the limit of 
large p. 

It is illuminating to consider the particular case where the V(x) is a screened 
Coulomb potential (Yukawa interaction) in two dimensions. We have 

V 2 V(x) - p 2 V(x) = -2tt5(x) (37) 

where is the screening length. Writing \x\ = r, V(x) = V(r) has the asymptotic 
behaviour 

V(r)«(-"; (r » .,_ r<< "" (38) 
[ exp(—pr)/y/JIr r » p 

From this we see that the integral 

d 2 x (cosh (PV(x)) - 1) (39) 



is convergent for large r when p > and for small r when (3 < 2. The importance of 
the screening is now revealed since its removal, by setting p = 0, leads to the integral 
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Bexpi — ^- (41) 



becoming infra-red divergent thus forcing the diffusion constant to vanish. It would be 
interesting to investigate this regime separately since the actual behaviour of the sample 
would presumably be sub-diffusive and hence scale dependent. 

The integral also diverges and the diffusivity goes to zero as (3 approaches 2 from 
below, that is as the temperature cools to T = 1/2. It is possible to analyse qualitative 
nature of the divergence since it is contained in the integral 

d 2 x (cosh((3V(x)) - 1) ~ / dr Ar 1 ' 13 ~ , (40) 

r<r JO 2 — (3 

for some A > and sufficiently small tq . It follows that the diffusivity behaves as 

4 9) o ( p a 

—2— pa B exp — 

which is of the form of the Volgel-Fulcher-Tammann law reported in fragile glass formers 
as they approach T g . It is important to recognise that the effect is one associated with 
short distances. A qualitative light is shed on the circumstances by considering the 
equilibrium tracer particle density in the presence of the scattering centres. At any time 
only the nearest such scatterer dominates the tracer particle probability density. If this 
has the opposite sign of charge to the tracer particle then the probability density is h(x) 
where 

h( \ = exp ((3V (x)) 
{X) f v d D x'exp{(3V{x')) ' 1 ' 

where we have placed the scatterer at the origin. Obviously this exhibits the same 
divergence as before and the density only exists for (3 < 2 . The physical interpretation 
of this result is that the tracer particle eventually becomes trapped at the point x = 
and thus diffusion is stopped. 

We now consider a system where the configuration of the scatterers is that obtained 
by freezing an equilibrium configuration of mobile interacting scatterers at equilibrium 
at inverse temperature (3, the same inverse temperature as the tracer particle i.e. with 
a distribution of the x n given by equation (JTTJ) with (3 = (3'. This situation would apply 
to a very mobile tracer diffusing in a background of very slowly moving charges of the 
same valence, for example an electron in a background of immobile ions. To facilitate 
the calculation of k we pass to the grand canonical ensemble for the scatters to obtain 

X 



N+,N- 



NANA V 



z n ++ n. exp g mV ( x . _ x .) J exp ^ qiV( Xi - x)j , (43) 

where the integration over x above is the spatial averaging of k(x) over V, 5 is the 
grand partition function and z is the fugacity of the positive and negative charges. The 
above expression can clearly be written as 
1 



K 



- £ (N„ + l)z N + +N - +1 Z N+ , N _ +1 , (44) 



N+,N. 
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where Z^ + ^_ is the canonical partition function for a system of N + /_ positive /negative 
charges. In the thermodynamic limit we thus obtain 
P_ 

2z ' 

which yields 



K 



(45) 



,(<?) 



Kq 



2z 
P 



For the case of a pure Coulomb interaction i. e. when \i 
V{x) = -ln(J 



\x\ 

T 



(46) 



(47) 



where L is a length scale which from here on we set to be 1. The statistical mechanics 
of the two dimensional Coulomb gas has been recently exactly solved ^Hj in the region 
where the model is thermodynamically stable (i.e. (3 < 2). For (3 > 2 the system 
is thermodynamically unstable and collapses unless a hard core interaction is included 
between the particles; this instability is present for the reasons mentioned previously. 
In the stable region it has been shown 



P 



1-/3/4 




0/4) 



r(i + /3/4) 

The resulting behaviour for k^ 9 ' when j3 <C 1 is 




1-/3/4 



.(48) 



Kq 



0>P) 



-/3/4 



exp 



2 7 + ln( 



7T 



4 



6 



C(3) 




C(3) 



where the the hyper-geometric function. Near the collapse point one finds 



/3 



K 



P7T 



,(49) 



(50) 



as /3 — > 2~. The interesting thing about these results is that the diflusivity is non- 
zero in all the region of stability, in contrast to the case of uniformly distributed 
Coulomb scatterers where it is zero in the thermodynamic limit in this region. Physically 
the fact that the Coulomb scatterers are allowed to equilibrate generates a screening 
characterised by a Debye length and the long range fluctuations of the electrostatic 
field are thus suppressed allowing a normal diffusion. The way in which the diffusion 
constant vanishes in equation (|5()j) as (3 — >• 2~ is also clearly very different to the case 
of a uniform distribution of Yukawa scatterers. 



6. Conclusion 



We have examined the problem of calculating the effective diflusivity of a Langevin 
particle diffusing in the force field generated by charged scatterers in two dimensions. 
The fact that the system is statistically electro-neutral means that the potential has 
the same distribution as —0. Using a general exact result connecting the effective 
diffusions constant of the gradient flow and varying permeability problem and an exact 
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result in two dimensions for the permeability problem allows us to solve exactly the 
problem of calculating k { 3) in the cases studied here. 

The results obtained give insight into the circumstances under which the random 
field may be taken to be Gaussian for the purposes of calculating the effective diffusion 
constant. It is clear from this study that the random field due to frozen scatterers will 
lead to transport properties that are generally different to that of a Gaussian random 
field with the same two point correlation function. It was also shown that the nature of 
the correlations between the point scatterers can drastically modify the behaviour of k[ 9 \ 
For instance the diffusion constant for a system electro-neutral uniformly distributed 
Coulomb scatterers is always zero at finite temperature, however if they are taken to 
in an equilibrium configuration of the same temperature as the tracer particles then 
4 9) is finite at high temperature. The exact results here should be useful in developing 
perturbative techniques to treat problems in higher dimensions and in two dimensions 
where there is only one type of scattering particle, i.e. no charge. 
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